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In a Comment [cond-m at/0506138 on our recent e-print [cond-m at/0505106 Liebsch claimed 
"excellent correspondence" between our high-precision quantum Monte-Carlo (QMC) data for the 
anisotropic two-band Hubbard model with Ising type exchange couplings and his earlier QMC 
results. Liebsch also claimed that the sequence of two orbital-selective Mott transitions, identified 
by us in this model, had already been reported in his earlier work. Here we demonstrate that both 
claims are incorrect. We establish that Liebsch's previous QMC estimates for the quasiparticle 
weight Z have relative errors exceeding 100% near transitions and cannot be used to infer the 
existence of a second Mott transition (for U C 2 ~ 2.5). We further show that Liebsch's attribution 
of our findings to his own earlier work is disproved by the published record. Consequently, the 
Comment is unwarranted; all results and formulations of our e-print remain valid. 



PACS numbers: 71.30.+h, 71.10.Fd, 71.27.+a 

In a recent e-print [l| we studied the anisotropic de- 
generate two-orbital Hubbard model (notation as in pj) 
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using high-precision quantum Monte Carlo (QMC) sim- 
ulations and showed (in accordance with independent 
slave-spin calculations |2() that this model contains two 
consecutive metal-insulator transitions and, hence, de- 
scribes an "orbital-selective Mott transition" (OSMT) 
0. This findingis to be contrasted with earlier results 
by Liebsch who reported a single Mott transition 

and explicitly excluded the occurrence of an OSMT for 
the same model. Our correction of this earlier picture 
established Q as a minimal model for the OSMT, seen 
experimentally in Ca2- a; Sr a ;Ru04, p| and demonstrated 
that the inclusion of spin-flip and pair-hopping terms (see 
and references therein) is not essential in this respect. 

While the high accuracy of our QMC data Q and the 
validity of our conclusions have not been challenged in 
any publication or e-print, Liebsch took exception to our 
discussion of his earlier work, first in || and then, in 
the form of a Comment on our work, in 9]. Liebsch 
now claims 0, that he not only observed the second 
transition but also reported on it in his earlier work 
and that, when he systematically mentioned only one 
single transition, he in fact meant two (one first-order 
and one second-order transition). Liebsch also claims |9j 
that the numerical data in our e-print hi are in "excellent 
correspondence" with his earlier work [6j , the implication 
being that any transition seen by us must necessarily also 
be contained in Liebsch's earlier work. 

The purpose of this Reply to Liebsch's Comment is 
twofold. First we show by a detailed comparison of both 
sets of QMC results that the data of are of insufficient 
quality near U C 2 and, hence, do not allow for any state- 



ment concerning the occurrence or non-occurrence of a 
second transition. In this crucial interaction regime near 
U C 2 there is, therefore, no "excellent correspondence" be- 
tween both data sets. Secondly, we demonstrate on the 
basis of explicit quotations from the published record 
that Liebsch did not identify a second transition in his 
earlier work. In fact he is generally viewed as the propo- 
nent of the single-transition scenario for this model and 
publicly attempted to adjust 0, Q this view only after 
submission of our paper [l| . In addition to these two main 
issues, we comment on the quality of Liebsch's QMC data 
near the first transition at U c \ and his quasiparticle spec- 
tra, and we detail our reanalysis ij of Liebsch's earlier 
IPT results @. 

Comparison of QMC data - In the following, we com- 
pare our QMC estimates for the orbital-dependent quasi- 
particle weights Z w idc, ^narrow Q with corresponding 
data by Liebsch [|| as well as with new numerically exact 
QMC data [n]]. We show that (i) our published QMC 
data is extremely accurate, with typical relative errors of 
1CP 2 , (ii) a careful analysis of this data reveals a second 
Mott transition at U C 2 ~ 2.5, and (iii) Liebsch's QMC 
data is much too noisy in the OSMT region - with rel- 
ative errors exceeding 100% near both transitions - to 
detect the second transition. In Fig.^,, the lower/upper 
sets of symbols and curves correspond to quasiparticle 
weights for the narrow/wide band. Our QMC results for 
discretizations At = 0.4 and At = 0.32 (circles/crosses) 
are practically on top of each other; a first kink in these 
estimates for both Z w i^ e and 2 narrow clearly signals a 
transition at U c \ « 2.0. Even on this scale (much larger 
than in 0, [tj), the significance of a change of slope 
in Z W ido in the region 2.4 < U < 2.7 is unclear. Only 
the analysis shown in Fig. ^> using a linear offset (as 
used for Figs. 1 and 5 in [l|) reveals a well- localized kink 
and, thus, a second transition at U C 2 ~ 2.5. Now con- 
sider Liebsch's QMC data (squares in Fig. already 
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FIG. 1: Comparison of discrete quasiparticle weights Z at 
T — 1/32: a) high- precision QMC data with minimal At 
dependence [jj (circles, crosses) clearly shows kinks in Z for 
narrow/wide band (lower/upper set of curves and symbols) 
at U c \ ~ 2.0; Liebsch's QMC data |(| (squares) deviates 
markedly near U c i and for U ~ 2.6. b) A subtle kink in 
the QMC data of Ref. Q] - clearly visible only after adding a 
linear term - indicates a second transition at U C 2 ~ 2.5; this 
signal is lost in the noise of the data of Ref. |6J . 

on the scale of Fig. QJi, it deviates strongly from our 
high-precision results near both transition regions. In 
particular, it cannot be used to detect the subtle transi- 
tion at U C 2, as demonstrated in Fig. ^d: all changes in 
slope of Liebsch's estimates for -Z w idc (squares and dashed 
line) are dominated by numerical errors; this data even 
displays an unphysical negative curvature of -Z w ide near 
U = 2.6. Liebsch's data also suffers from too coarse a 
grid and, in particular, from the lack of data points be- 
yond the retroactively claimed [|| transition. We further 
stress that Liebsch's criterion that iv w ide "becomes very 
small near Ub ~ 2.7 eV" Q is not a valid criterion for a 
transition. 

In order to fully quantify the numerical errors (includ- 
ing Monte Carlo statistical errors, convergence errors and 
discretization errors) of all data sets shown above, we 
have extended our QMC simulations to smaller dis- 
cretizations (At = 0.25 and At = 0.20) and have ob- 
tained numerically exact reference data by explicit extra- 
polations At — > 0. Relative errors computed by compar- 
ison with this reference data are shown in Fig.|3for both 
the wide and narrow band (lower and upper panel, re- 
spectively). The relative errors in our QMC data (circles, 
crosses) are visible only at the narrow-band Mott transi- 
tion (in ^narrow for U = 2.0); throughout the rest of the 
OSMT region, our relative errors are typically smaller 
than 10 -2 . The errors of Liebsch's data (squares in Fig. 
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FIG. 2: Estimates of relative errors in Z at T = 1/32: Lieb- 
sch's QMC results (squares) deviate significantly from new 
numerically exact QMC data - bymore than 100% near tran- 
sitions. In contrast, already our [l| finite-Ar QMC data (cir- 
cles, crosses) is accurate within about 1% (except for Z narrow 
at U = 2.0). 

|2J) are in this region two orders of magnitude larger and 
exceed 100% at both transitions (as identified by us); 
only in the metallic phase for U < 1.6 (not included in 
this comparison), Liebsch's data appear reasonably accu- 
rate. We conclude that the second transition, identified 
by us, cannot be seen in Liebsch's QMC estimates of Z. 
This already disqualifies the Comment. 

Liebsch's published record on OSMT - We now come 
to the second part of our Reply, in which we show on 
the basis of the published record, that Liebsch in 
fact proposed a scenario very different from the OSMT 
scenario found by us, namely the existence of only a sin- 
gle Mott transition. Indeed , up to the recent surge of 
e-prints on the subject, [ll 0, HH u2, LLil Liebsch was gen- 
erally recognized as the proponent of the single-transition 
scenario, as witnessed, e. g., b y various recent statements 
to this effect in Refs. Hill E Q. It was only after 
submission of our paper that Liebsch 0, |j| attempted to 
reinterpret his earlier work 0, 0, @ ■ That such a rein- 
terpretation is in fact precluded by the published record 
will now be shown. 

Chronologically Liebsch's first paper on the anisotropic 
degenerate multi-orbital Hubbard model is , which car- 
ried the title "Absence of orbital-dependent Mott transi- 
tion in C '&i - x Sr x RuO . The abstract informs us, that 
the results suggest "a common metal-insulator transition 
for all t2 g bands at the same critical IP'. Liebsch ex- 
plains on the basis of general arguments from the theory 
of phase transitions (his Ref. 17) that in several "con- 
ceptually closely related" problems "a true separation 
of phases does not seem possible" . In accordance with 
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this expectation, Liebsch then finds, using QMC/DMFT, 
"a behavior that differs fundamentally from the one re- 
ported" by Anisimov et al. [3j, characterized by "one 
common metal-insulator transition [existing] at an inter- 
mediate critical value of U between those of the isolated 
d X z,yz and d xy bands". Needless to say that this finding 
is irreconcilable with the existence of two separate tran- 
sitions. In 0, Liebsch calculates the density of states of 
Si'2Ru04 for various parameter values and temperatures 
and deduces from the results the existence of "a common 
Mott transition" satisfying a U^ yz < < U^f, i.e., 

the Mott transition for the actual ti g complex requires 
a critical U between those of the independent d xz , yz and 
d xy bands" . Again, this leaves no room for two separate 
Mott transitions. 

Chronologically next in line is Ref. jjj, the abstract 
of which informs us, that "interorbital Coulomb interac- 
tions in nonisotropic multiorbital materials give rise to a 
single Mott transition" , implying (according to the third 
paragraph of the text) "all subbands to be either metal- 
lic or insulating" . The text then goes on: "The critical 
Coulomb energy U c of the interacting system lies between 
the U C i of the isolated subbands. Thus, narrow subbands 
are less correlated and wide subbands are more correlated 
than in the absence of interorbital Coulomb interactions. 
Coexistence of metallic and insulating behavior in differ- 
ent subbands does not occur." Similar statements are 
repeated both in the bulk of the text and in the captions 
to Figs. 1 and 2 of that paper, leading to the conclusion 
that "the present picture does not support the 'orbital- 
selective' Mott transitions". 

Last in line of the early works of Liebsch on this subject 
is Ref. 0, which has as its title "Single Mott transition 
in the multiorbital Hubbard model", while the abstract 
states that both the IPT approximation and QMC "give 
rise to a single first-order transition rather than a se- 
quence of orbital-selective transitions" . This conclusion 
of Ref. 0, in combination with the alleged "complex" 
finite temperature behavior of the subbands, constitutes 
a serious misinterpretation of the orbital-selective Mott 
transition, which actually occurs in this model. Ref. 6] 
states (in Section II) that the quasiparticle weights ob- 
tained from IPT "exhibit first-order transitions at pre- 
cisely the same critical U C (T). This picture supports our 
previous results suggesting a single Mott transition in a 
nonisotropic multiorbital environment, in contrast to the 
orbital selective transitions found [by Anisimov et al. and 
Koga et al.]." We shall see below that in fact IPT, too, 
predicts an orbital-selective Mott transition. Using the 
QMC-DMFT, Liebsch then shows (in Section III) that, 
"as in the case of IPT, there is no evidence for a second 
transition in the wide subband at larger U" . The paper 
concludes (in Section IV), that "the present results con- 
firm our previous finding, namely, the existence of a sin- 
gle transition rather than a sequence of orbital-selective 
transitions as the on-site Coulomb energy is increased" . 



Thus the occurrence of an OSMT is explicitly denied in 
both Ref. 5] and Ref. @. 

We now consider Liebsch's Comment [t| . In the Com- 
ment, Liebsch ignores his previous work 0, Q and fo- 
cuses on his QMC calculations at T = 0.031 eV in [||. 
This focus on Ref. Q is remarkable, since Ref. 0, as 
we saw above, just confirms his "previous finding" . In 
it is then stated that in Q two transition regions (at 
U a and Ub) were identified, with bad-metal behavior and 
a breakdown of Fermi-liquid behavior for the wide band 
in the "intermediate phase" U a < U < Ub- Here we 
wish to record that in Q neither the notation U a nor 
Ub was used, that "Fermi-liquid behavior" (or the break- 
down thereof) was not mentioned in the entire paper, 
that Liebsch in fact argued (by comparison with data for 
T = 0.02 eV) that "both subbands undergo a common 
transition at the same Coulomb energy and that the wide 
subband exhibits pronounced bad-metal behavior above 
U c " and that he concluded (see also above) that "there is 
no evidence for a second transition in the wide subband 
at larger W . We further note that, since there is no "sec- 
ond transition" and no "J7 b " in [||, the hysteresis behav- 
ior cannot have suggested that "only the lower transition 
is first-order" . In Q Liebsch states that he argued in Q 
"that the Mott transition in the nonisotropic Hubbard 
model is governed by a single first-order transition . . . " , 
suggesting that he left room for a second continuous tran- 
sition at some larger value of U . Comparison with the 
above quotations from the published record shows that 
this interpretation of Ref. |g] is untenable. 

We briefly comment on Ref. .8] , in which our work Q is 
also cited. While ; g] is, generally speaking, an interesting 
extension of earlier work [Tr| to finite temperatures, this 
preprint, too, contains several inaccuracies. For instance, 
it is not true that (as claimed in the second column) 
the phase diagram of || is in agreement with that of 
@ (cf. Fig. 5 in and Fig. 4 in @), or that (fourth 
column) the two-transition scenario of Fig. 2 "confirms 
the picture" of Ref. Q , or that (sixth column) Liebsch's 
use of "single Mott transition" [|| did not "imply the 
non-existence of the continuous transition . . . " . All of 
this is irreconcilable with the passages from the published 
record, quoted above. We further note that "breakdown 
of Fermi- liquid behavior" occurs in this paper for the first 
time, not in Q (as claimed in 

Comment on Liebsch 's spectra - Quasiparticle spectra 
play a subordinate role in Liebsch's work 0,13,0], accord- 
ing to Ref. in order to "avoid the uncertainties from 
the maximum entropy reconstruction" . Here we discuss 
them nonetheless, primarily because our spectra Q] are 
compared to those of Ref. 0] in 0. Generally, Lieb- 
sch's spectra (obtained from QMC using the maximum 
entropy method) appear overly broad with strong high- 
frequency tails and very smooth behavior at low frequen- 
cies. This lack of resolution obscures the observation of 
gaps: e.g., the narrow-band spectrum for U = 2.1 (solid 
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FIG. 3: Reanalysis of IPT estimates for Z: A high-resolution 
plot of Zwide (extracted from the preprint to jjjj ) reveals tran- 
sitions near U C 2 ~ 3.6 — 3.7 as kinks in otherwise linear 
fits which represent the data very well for T < 0.02. The 
rounding-off at T = 0.03 (and U = 2.6) - if genuine - may 
indicate that Tip T < 0.03. Insets: the transitions apparent in 
the main panel are obscured in representations similar to 0,0 
- on this scale, the IPT data is consistent with featureless fits 
(e.g., of the form aexp[— bll c ] - dashed lines). 

line in upper- most panel of Fig. 2 in $j) appears as still 
metallic [with N{lu = 0) « 0.022, more than 3% of the 
noninteracting value] ; a true gap for this interaction [with 
N(u) = Q) = 0] has only been established in our work (cf. 
Fig. 4 in 0). In Liebsch's data, this gap in the narrow 
band is only well-established at U = 2.4 (second panel of 
Fig. 2 in [9j). The lack of a visible gap in the wide-band 
spectrum is attributed in to the filling of a pre-existing 
gap at finite temperatures: "Only if U is increased to 
about 2.7 eV is the gap in the wide band large enough to 
not be obliterated by this temperature broadening." The 
delayed opening of gaps (with respect to the "single Mott 
transition" pinpointed at U — 2.1) is explicitly treated 
on equal footing for both bands: "For U = 2.7 cV ... a 
gap opens up [in the wide band]. The narrow band un- 
dergoes a similar cross-over behavior, except at slightly 
lower values of U." @ We conclude that the broadness 
of Liebsch's spectra (which must be attributed to noisy 
QMC data) prevents a proper characterization of phases 
as metallic or insulating; in addition, more grid points 
would have been necessary in order to see the second 
transition. In comparison, our spectra show well-defined 
gaps with sharp edges (Fig. 3 in [J and lower three panels 
of Fig. 2 in 0); still, we would not have relied on spectra 
alone for the identification of the second transition. 

Reanalysis of IPT-data of Ref. 0/ - Finally we prove 
our previous statement that, in contrast to Liebsch's find- 
ings in 0.0]. the IPT approximation, too, describes two 
separate transitions rather than "a single Mott transi- 



tion" 6]. For this purpose we use Liebsch's own IPT- 
data [6j for the quasi-particle weight, which is plotted 
in Fig. The insets show data for T = 0.02 and 
T = 0.03 for both the narrow and the wide band. The 
main panel shows Z w id e (U) for various temperatures and 
3.2 < U < 4.0. As can be seen from Fig. [31 the be- 
havior of Z w id c (U) is approximately piecewise linear on 
both sides of U C 2(T) ~ 3.6 — 3.7, where U C 2(T) increases 
slightly with decreasing temperature. This second tran- 
sition at U C 2 is not visible in the small-scale graphical 
representation of the IPT data in Figs. 1-4 of [7] (similar 
to the insets in Fig. [3| and was mentioned by Liebsch 
only after he was informed of its existence by us . 

Summary - We have shown that, in contrast to claims 
in H, there is no "excellent correspondence" between 
the numerical data of our paper pj and those of Ref. 0] , 
and that the conclusions reached in Q indeed correct 
previous erroneous statements concerning "single Mott 
transitions" in Q. Thepublished record clearly demon- 
strates that Liebsch |4j,|5j,|g] did not observe two separate 
transitions. A comparison of Liebsch's QMC data with 
our high-precision results further shows that the second 
transition was hidden in by huge error bars (of up to 
100% in the critical region) and inadequate data analysis. 

We thank E. Jeckelmann and D. Vollhardt for discus- 
sions and acknowledge support by the NIC Julich and by 
the DFG (Forschergruppe 559, B1775/1). 
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